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ABSTRACT 

By comparing the outcome of N-hody calculations that include primordial residual- 
gas expulsion with the observed properties of 20 Galactic globular clusters (GCs) for 
which the stellar mass function (MF) has been measured, we constrain the time-scale 
over which the gas of their embedded cluster counterparts must have been removed, 
the star formation efficiency the progenitor cloud must have had and the strength of 
the tidal-field the clusters must have formed in. The three parameters determine the 
expansion and mass-loss during residual-gas expulsion. After applying corrections for 
stellar and dynamical evolution we find birth cluster masses, sizes and densities for the 
GC sample and the same quantities for the progenitor gas clouds. The pre-cluster cloud 
core masses were between 10 5 — 1O 7 M0 and half-mass radii were typically below 1 pc 
and reach down to 0.2 pc. We show that the low-mass present day MF (PDMF) slope, 
initial half-mass radius and initial density of clusters correlates with cluster metallicity, 
unmasking metallicity as an important parameter driving cluster formation and the 
gas expulsion process. This work predicts that PD low-concentration clusters should 
have a higher binary fraction than PD high-concentration clusters. 

Since the oldest GCs are early residuals from the formation of the Milky Way 
(MW) and the derived initial conditions probe the environment in which the clusters 
formed, we use the results as a new tool to study the formation of the inner GC 
system of the Galaxy. We achieve time-resolved insight into the evolution of the pre- 
MW gas cloud on short time-scales (a few hundred Myr) via cluster metallicities. The 
results are shown to be consistent with a contracting and self-gravitating cloud in 
which fluctuations in the pre-MW potential grow with time. An initially relatively 
smooth tidal-field evolved into a grainy potential within a dynamical time-scale of the 
collapsing cloud. 
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1 INTRODUCTION 

Globular clusters (GCs) and old low-mass stars have often 
been used as local probes of Galaxy formation, since they 
preserve information about ancient times. In the past, the 
formation of the Milky Way (MW) has been investigated by 
means of kinematic studies of stars and their abundances 
dCarorio et al.ll2007l ; iBeers et al.ll2002l ; iBeers fc Chiball200ll ; 
lEggen et all 1 19621) . as well as by horizontal branch mor- 
phology, metallicity and kinematic measurements of star 
clusters dBekki et al.ll2007l ; IMackev fc van den Berghll2005l ; 
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IMackev fc Gilmord 12004 IZinnl Il993l ; ISearle fc Zinnl Il978l ). 
In terms of the origin of the GCs, this led to a picture of 
Galactic formation in which GCs are divided in the old halo 
(OH) and young halo (YH) clusters. The OH clusters are lo- 
cated inside a Galactocentric distance of doc ~ 8 — 10 kpc. 
Many of them appear to have formed coevally with the col- 
lapse of the protogalaxy (lEggen~et al. 1962; Salaris & Weiss! 
l2002l ; lDe Angeli et al.ll2005l ; iMarfnTYanch et ai1l2009T l. The 
YH clusters have rfGC 8 — 10 kpc an d have been accreted 
over several Gyr (|Searle fc Zinnlll978l ). All this information 
was gained from knowledge about the present day (PD) pa- 
rameters of GCs. Knowing about the initial conditions at 
star cluster birth, however, would provide a deeper insight 
into the early formation processes of the MW since it would 
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allow to probe directly the e nvironment in which the GCs 
formed l|Goodwinlll997l Il998l ). 

Cluster masses at birth were larger than they are today 
since clusters suffer mass loss due to primordial residual- 
gas expulsion, stellar and dynamical evolution. And since 
the expulsion of gas leads to subsequent expansion, young 
and gas-embedded clusters were much more compact and 
denser than they are nowadays. Star clusters are the PD 
gravitationally bound remnants of these dense o bjects af- 
ter the stars emerged from their natal cloud (iTutukov 
197S ; iKroupa. Aarseth fc Hurlevll200ll ; IBastian fc Goodwin 
20061 : iBaumgardt fe Kroupal l2007f ). The overall change in 
the potential due to the gas loss leads to cluster expan- 
sion and the loss of stars, with the initial conditions at 
the onset of this process deciding about cluster survival 
or destruction. The majority of the freshly hatched clus- 
ters are destroyed during this violent phase dLada fc Ladal 
20031; IBastian fc Goodwin] 120061; Ide Griis fc Goodwin! 1200a 
Bastianl 120081 ; iGieles fc Bastianl [2008) and their member 
stars become stars of the field. In this view, the OH GCs 
are the massive remnants of an initial population of embed- 
ded clust ers that rapidly formed the population II halo of 
the MW (IKroupa. Aarseth fc Hurlevll200ll ; IKroupa fc Bolivl 



120021 ; IBaumgardt. Kroupa fc Parmentien l2008l ). Some of 
them ended as bound clusters which, after stellar and two- 
body relaxation driven dynamical evolution, we can observe 
nowadays. So in order to understand physical properties of 
star clusters today it is essential to understand how star clus- 
ters formed, because the birth configuration determines the 
fate of a star cluster. 

The time, tm, over which the natal gas is 
removed from the cluster determines crucially 
whether a star cluster survives gas expulsion or not. 
IBaumgardt. Kroupa fc Parmentierl 
analytic formula to calculate tm, 



(2008) provide an 
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based on the amount of energy needed to be put into the 
gas to overcome its potential energy. The deeper the poten- 
tial, i.e. the larger the progenitor cloud mass, M c \, the more 
difficult it is to remove the gas. A large star formation effi- 
ciency (SFE), e, leads to more and also more massive stars 
with stronger winds and radiation and there is less gas to re- 
move. Finally, the larger the half-mass radius, th, the faster 
the gas is expelled since the overall potential is shallower for 
fixed Mci and e. 

The gas expulsion time, thus, depends on the mass and 
the size of the cluster. From theoretical considerations mass 
and radius are related. Dependent on the exact form of the 
mass-radius relation of young star clusters, the gas expulsion 
time-scale, tm, is affected more strongly or weakly. Virialised 
gas cores (radius r c and core mass m c ) are expected and ob- 
served to show a strong mass-radius relati on, scaling as r e oc 
jnV 2 (i-e. constant surface density, e.g. lHarris fc Pudritj 
1994). Therefore the mass-radius relation of young clusters 
is expected to display a similar behaviour. If this relation is 
valid the gas removal time-s cales would be essentially m ass 
independent, tm oc mj (|Parmentier fc Fritz 3 120091 '). In 
contrast, observations of young clusters do not sh ow any 
signifi cant mass-rad ius relation, r c ; oc rn®~ ' 1 (e.g. iLarsenl 
120041 ; iKroupal 120051 ). and the influence on the gas expul- 
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Figure 1. Concentration parameter, c = log 10 (rt/r c ), vs. low- 
mass MF slope, a, (the DMPP plane). Weakly concentrated clus- 
ters are strongly depleted in low-mass stars and no cluster with 
a high concentration and a depleted MF is found (filled dots). 
This trend (black solid line, eq. [3J| can't be understood in terms 
of purely secular dynamical evolution. However, JV-body integra- 
tions of mass-segregated clusters at the time of the emergence 
from their birth molecular cloud with unresolved binaries in them 
(squares with error bars) reasonably reproduce the observed trend 
within the observational limits (dashed lines). 



sion time-scale is a strong er function of mass, tm oc m}J 2 
|Parmentier fc Fritzdl2009r ). 

Obs ervationally determ ined SFEs lie between e = 0.2 
and 0.4 (jLada fc Ladall2003l ). so typically most of the mass 
remains in the gas to be expelled. JV-body modelling showed 
that SFEs can be as low as 10 — 20 per cent if the gas 
expulsion time-scales are sufficiently long and has to be 
at least 33 per cent if the g as is lost instantaneously in 
order for a cluster to su r vive (|Lada. Margulis fc Dearborn! 
1984 iGever fc Burkertl I2001I ; iBoilv fc Kroupal l2003al lbl; 
Baumgardt fc Kroupal 120071 ). SFEs may vary locally and 
could be different between high- and low-mass cores, which 
is also a viable explanation t o wipe out the mass-ra dius re- 
lations observed in gas cores (|Ashman fc Zepdl200ll ). 

In Sec. [2] we present the observational data to be com- 
pared to the results of JV-body experiments and in Sec. [3] we 
explain how we model the initial cluster size, mass and den- 
sity. We derive the initial conditions and discuss correlations 
between and among initial and present day cluster parame- 
ters in Sec. [4] In Sec.[5]we connect the primordial conditions 
constrained in the former sections to develop a picture for 
the assembly of the old population of GCs located in the 
inner halo. Sec. [6] summarises our main results. 
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Figure 2. Global cluster metallicity, [Fe/H] dHarrid[l996f) , vs. 
low-mass MF slope, a, as in Fig. [l] Data points are labelled with 
the respective NGC or Pal catalogue number of the GC. The the- 
oretical expectation would be that otherwise similar systems with 
a lower metallicity ought to have a bottom-light IMF compared to 
more metal-rich systems. Instead, clusters having a larger metal- 
licity ([Fe/H] > —1.5) are depleted in low-mass stars (a < 1) 
giving support to the gas expulsion scenario and arguing against 
a variable IMF (see the text). 



2 THE DE MARCHI DIAGRAM AND THE 
METALLICITY-MF SLOPE RELATION 

The number of stars in the mass interval m, m + dm is 
dN — £(m)dm, where £(m) oc m~ a is the stellar mass- 
functi on (MF) with index a. |Pe Marchi. Paresce fc Pulond 
|2007l . hereafter DMPP) showed that for Galactic GCs the 
value of a in the mass range 0.3—0.8 Mq correlates with their 
concentration parameter c = log (r t /r c ), i.e. the logarithmic 
ratio of tidal- and core-radius (Fig. [TJ filled circles). The be- 
haviour of a with c cannot be understood to be the result 
of secular dynamical e volution without further assumptions 



mpt 

(for details see DMPP;lMarks. Kroupa fc Baumgardtlhoosl ; 
iBaumgardt, De Marchi fc Kroupal [20081 ) . By fitting by eye 
DMPP proposed that the clusters obey a relation of the form 



q(c) 



2.3 



2.5 



(2) 



(solid black line in Fig. [T] hereafter referred to as the DMPP 
relation). 

One exciting explanation for this trend is that it may 
signify for the first time evidence for a variability of the ini- 
tial MF (IMF) at low masses. Low-c clusters would form 
fewer low-mass stars than highly concentrated clusters. One 
would naturally expect this to depend on the metallicity, 
[Fe/H], of the cluster. Fig. f2] displays the observed global 
MF slope versus the metallicity of the clusters in the sample. 
The anti-correlation between the two quantities is signifi- 
cant (see Tab. [2] for a compilation of Spearman rank order 
correlation coefficients) and appears even slightly stronger 
when neglecting the four young, high- metallicity disc clus- 



ters NGC 104, 6352, 6496 and 6838. The data appear to 
suggest that stellar IMFs need to be flatter (smaller a) in 
nowadays metal-rich clusters than in their metal-poor coun- 
terparts. In low metallicity environments, however, star for- 
mation should be terminated later by radiative feedback be- 
cause the radiation couples less effectively to metal-poor gas. 
Stars would accrete matter over a longer time such that low- 
mass stars should be under-abundant (|Adams fc Fatuzzol 
1996). Additionally, the Jeans mass is expected to be larger 
because of less effective cooling in gas clouds, ag ain sug- 
gestin g the IMF to have a larger average mass (e.g. iLarsonl 
I1998T ). From these simple arguments we would expect to 
find depleted MFs in low- metallicity clusters, which is not 
observed. For these reasons the data in Fig. [J| are incom- 
patible with our theoretical knowledge of star formation and 
up to now no purel y star-formation-bas ed explanati on for 
this trend exists (seelMcClure et al Jl 198(1: ISmith fc McClurd 
ll987l ; lDiorgovski et al.lll993D . 

In an investigatio n of these issues, 

iMarks. Kroupa fc Baumgardtl (|2008l . hereafter MKB08) 
showed, usin g TV-body integrations of you ng, gas embedded 
clusters from lBaumgardt fc Kroupal (|2007| . hereafter BK07), 
that the loss of the primordial residual-gas from initially 
mass-segreg ated clusters starting with a canonical IMF 
|Kroupal l200lfl and containing unresolved binaries (Fig. [T] 
squares with error bars) reasonably reproduces the trend 
within the observational limits (dashed lines). 

Wtthin the framework of gas expulsion we can now 
qualitatively explain the metallicity trend seen in Fig. [J| 
Stellar winds and the radiation of stars will lead to de- 
struction of the gas cloud, in which the newly formed 
stars are embedded. Just as for the metal-dependend mass 
loss rates in stars (iMokiem et~aT1 120071 ; IVink et all 1200 ll ; 
iKudritzki fc Pulsll2000l ). we suggest that the radiation of 
stars, which formed from a high-metallicity cloud, couples 
better to the gas through the dust, so that the residual- 
gas can be expelled more quickly via the efficient deposition 
of momentum into the interstellar medium. This leads to 
faster gas blow-out in high-metallicity clusters. In turn, the 
stronger change of the potential in a given time results in 
stronger expansion, i.e. a lower final concentration after gas 
expulsion, and stronger mass-loss over the tidal boundary, 
i.e. a lower value of a. This is observed by DMPP. 

However, cooli ng is more efficient in higher met allic- 
ity environments (|Larsonlll998l ; lAdams fc Fatuzzol 19961 ) and 
the emanating photon flux from the st ars is higher if t heir 
atmospheres have a lower metallicity (Kudritzki 20021 ). It 
is thus the interplay between the above three effects that 
determine the influence of metallicity on the gas expulsion 
process. We will demonstrate that metallicity affects the 
time-scale of gas expulsion (Fig. [B] below) so that our idea 
in favour of a [.Fe/-ff]-dependend gas loss process appears 
justified. 

This supports the picture that residual-gas loss is indeed 
the driving mechanism to initiate the relation in Fig. Q] and 
there may be no need to invoke a variable IMF at low masses. 



1 The canonical IMF has a 
a = 2.3 for m/ Mq > 0.5. 



1.3 for 0.08 s; m/M ^ 0.5 and 
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3 THE MODELS 

Our aim is to calculate the mass of the star forming cloud, 
M c i, its size, rn,a and average density within the half-mass 
radius, p a v,ci, each of the DMPP clusters must have formed 
from. This will enable us to investigate how GCs formed and 
what they and their birth sites looked like before residual-gas 
expulsion. 

For each GC we constrain a gas expulsion time-scale, 
tm /icross, in units of the cluster initial crossing time, a 
star formation efficiency (SFE), e, and a value for the ini- 
tial ratio of half-mass to tidal-radius, rh,ci/rt,ci, in Sec. 
|4] These are the model parameters in BK07 which de- 
termine the fraction of final to initial half-mass radius, 
/r(res.gas) = rh, t R n /rh,ecl, and the bound stellar mass- 
fraction, fja (res.gas) = Ma n /M ec i, after the gas loss process 
(their figs. 1 and 4). Here, M ec \ = eM c \ is the mass in stars 
in the embedded cluster that formed out of the gas cloud 
and rh,eci is the corresponding half-mass radius. We find a 
value for f r (res.gas) and ju (res.gas) for each of the 20 GCs 
in the DMPP sample by interpolating between the BK07 
model grid. 

We consider the PD and M as the result of three 
independent evolutionary steps. The PD half-mass radius, 



= /(dyn.evo.) x / (st. mass loss) , 



(9) 



?Xpd = /r(dyn. evo.) x / r (st. mass loss) x 
f r (res.gas) x r hiCC i , 



(3) 
(4) 



comes about from the change of the initial half-mass radius 
due to residual-gas expulsion, /,. (res.gas), from mass loss 
due to stellar evolution, /r-(st. mass loss), and due to the 
dynamical evolution over a Hubble-time, / r (dyn. evo.). We 
assume the latter to be negligible since the half-mass ra- 
dius is approximately constant during dynamical evolution 
ijKupper. Kroupa fc Baumgardtll2o"ol ). so / r (dyn. evo.) w 1 
and eq. Q reduces to 



f/i.pd — /r(st. mass loss) x / r (res.gas) x r^., 



(5) 



For these processes f r > 1 holds. According to 
iBaumgardt fc Makind (2003), clusters having a canonical 
IMF with stellar masses between 0.1 and 15 Mq initially 
loose about 30 per cent of their mass due to stellar evolu- 
tion of the massive stars. In clusters with a fully populated 
IMF (up to 150 Mq) this mass loss will be larger and closer 
to 40 per cent. The loss of this gas from the cluster gives 
rise to further cluster expansion, whic h we assum e to be 
adiabatic. In the case of slow expansion l|Hillslll980r i the PD 
half-mass radius is then given by 



rh,PD = ~~~~~~ => fr (st. mass loss) 
0.6 



?"ft,PD 



1 

06 ' 



(6) 



where Th.rin = f r (res.gas) x rh iCC i is the previously denned 
cluster half-mass radius after residual-gas expulsion and the 
factor 1/0.6 accounts for the mass loss. 

We do the same to constrain M ec j and we include the ef- 
fect of mass loss due to dynamical evolution. The PD stellar 
mass of a cluster is then 

Mpd = /m (dyn.evo.) x /m (st. mass loss) x (7) 
f M (res.gas) x M ecJ , (8) 

where /m (res.gas) is taken from the BK07 results as de- 
scribed above and /m < 1 for all processes. The product of 
the first two factors is 



Mpp 

M fln 

where Mfl n = /m (res.gas) x M cc \ is the above defined mass 
after residual-gas loss when the cl uster has re-virialised. Fo r 
Mfi n we use the values derived in IBaumgardt et all l|2008h . 
who calculated initial cluster masses for the DMPP sample 
considering dynamical and stellar evolution only. 

Via eqs. ((5} and |[8j we can calculate the initial half- 
mass radius, rh, cc i and stellar mass, M ec i for the young, gas 
embedded counterparts of the clusters in the DMPP sample. 
We then define the initial stellar density as 



Pav,ecl — Pav,ecl ( 



rh,ec\) 



Meci (< Th,cc\) 
V(rh,ecl) 



(10) 



i.e. the total stellar mass within rh,ec\ divided by the volume 
of the sphere with radius rh, cc \. 

The corresponding quantities for the star forming cloud 
are obtained by replacing M cc \ with M c \ and assuming 
Th,c\ = Vh,eci, i.e. stars and gas followed the same radial pro- 
file. From here on we set r^ci = ffe.ecl- Note that the above 
initial densities refer to the point of time when the forming 
cluster is close to virial equilibrium before gas expulsion and 
approximately corresponds to the state of maximum density 
after cluster cloud-core contraction. 



4 INITIAL CONDITIONS 

In order to derive the initial conditions described in Sec. [3] 
we first need to determine the time-scale of gas expulsion, 
tm /Across, the SFE, e, and the strength of the tidal- field, 
Th.cifrt.ci- In Figs.[3]and|4]we compare the results of MKB08 
with the observed data to constrain the initial rh, c i/r ttC i and 
tm /tcross values. MKB08 considered the effects of gas expul- 
sion and unresolved binaries but they didn't include stellar 
and dynamical evolution in their analysis. 

However, dynamical and stellar evolution do change a 
and c and we will discuss their influence later (Sec. 14. 5[) . 
At this point we only note that the core-radius, r c , en- 
tering the concentration parameter, c, is probably changed 
strongly by processes such as core collapse. The half-mass 
radiu s, r^, is much mor e stable against dynamical evolu- 
tion l|Kiipper et al.| [20C)8). In Figs. [3] to [5] we therefore use a 
modified concentration parameter, 



Cmod 



logi 



n 



(11) 



310 \rh P J 

and use it instead of c as a more robust measure. The ra- 
dius r^p = ?Xpd is the PD projected half-mass ra dius and 
r t is the PD tidal radius, both from iHarrisl (| 19961 . revision 
of 2003). The tidal radius is assumed not to vary strongly 
after the cluster re-virialised after gas expulsion. The trend 
between a and c mo d is similar to that in Fig. [T] and enables 
us to gain more reasonable results. 



4.1 Tidal- field strength, time-scale of gas 
expulsion and star formation efficiency 

From Fig. [3] we find that the primordial tidal-field grows as 
Cmod decreases. We grouped the initial tidal-field strength 
used in the models of BK07 into three groups. We draw 
boxes around areas containing (most of) the models with 
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Figure 3. As in Fig. [T] but with the modified concentration pa- 
rameter (eq. |1H on the abscissa. The theoretical data are coded 
for the initial tidal- field strength, r h c \/r t c y in the simulations 
of BK07 with different symbols. Open symbols and crosses depict 
the TV-body results, the filled circles labelled with their respective 
catalogue number are the data from the DMPP globular clus- 
ter sample. Quadrangles and triangles show areas with the same 
Th cl/ r t cl- AH initial models star ted from the dashed line, which 
is the slope of the canonical IMF (Kroupa 2001) in the mass-range 
0.3 - 0.8 M (a « 1.72), in which the MF slope of the observed 
GCs was measured. 



the same rh. c i/rt. c i-values by eye. Comparing the results of 
the integrations (open symbols and crosses) with the ob- 
servational data (filled circles) we estimate the strength of 
the tidal-field for each DMPP cluster. If we can't assign 
a unique number to a GC, i.e. if a cluster is located in 
an overlap region, we choose a maximum and a minimum 
tidal-field strength defined by the overlapping regions and 
assign the mean value. NGC 6752, for instance, is assigned 
'Xci/Vt.ci = 0.02, while for NGC 6656 r hyC \/r t , c i = 0.05 was 
chosen, i.e. the average of 0.02 (minimum tidal-field) and 
0.08 (maximum tidal-field). For GCs outside any of the re- 
gions (NGC 2298 and NGC 7078) the minimum and maxi- 
mum tidal-field is chosen according to the r^ci/r^ci-values 
occurring around the same value of c mo d- 

In Fig. [4] we divide the models into those with slow 
(rAf/tcross » 1), intermediate (rAf/tcross ~ 1) and fast 
(tm /tcross -C 1) gas removal. The location of a model clus- 
ter in the a — c plane is well described by a function of 
the form of the DMPP-relation, eq. [2l for slow and inter- 
mediate gas expulsion times. The trend is flat for fast gas 
removal. Areas in Fig. [4] containing models with similar gas 
expulsion times are enclosed by two lines of the same type. 
As for the tidal-fields we assign values for tm /tcross to each 
cluster in the sample corresponding to their position in the 
DMPP diagram. E.g., NGC 6809 gets r M /t CIOSS = 10, NGC 
5272 is assigned tm Across = 0.5 (~the average of 0.1 and 
Wcross) and, according to the models, NGC 6352 needed 5 
crossing-times to remove its residual-gas. 



Figure 4. Same figure and symbols as in Fig. [3] but now coded 
for the residual-gas expulsion time-scale, tm /tcross ■ Two lines of 
the same type enclose regions with the same tm / tcross • Adiabatic 
(TM Across 3> 1, dashed lines) and intermediate (jm /tcross ~ 1, 
solid lines) residual-gas expulsion times approximately follow the 
form of the DMPP relation, eq. [2] The trend for explosive gas 
removal (tm /tcross <C 1, short-dashed lines) is flat. 
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Figure 5. Same picture and symbols as in Figs. [3] and [4] but here 
coded for the SFE, e. SFEs of every value are distributed over the 
whole range of concentration except that for the depleted models 
( a ^ 1) mainly large SFEs (> 33 per cent) are found. High SFEs 
for low-a clusters are needed to retain a bound cluster. 
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The SFE shows no clear behaviour in the DMPP plane 
although it is an important parameter deciding about clus- 
ter survival or cluster destruction (Fig [5]). For the models 
that are still closely resembling their IMF (a > 1) basically 
all values from the SFE-range covered in the iV-body inte- 
grations occur. Solely the strongly depleted (a < 1) model 
clusters nearly show large SFEs only (e > 33 per cent). For 
these reasons it is difficult to make a reasonable guess for 
the SFE for any of the DMPP clusters. We decided not to 
directly compare the observations with the models as before, 
but to look which SFE is at least necessary for a given pa- 
rameter combination (fTt.clA't.cl, TM /Across) to retain a bound 
cluster. The SFE is then chos en such that at least four clus - 
ters from the grid of models of lBaumgardt fc Kroupal (|200?t ) 
around (r h ,ci/r t ,ci, tm Across) survive gas expulsion^. Hence, 
we find lower limits to the SFE, ei ow . 

4.2 Size, mass and density 

For our modelling in Sec. [3] we can now directly get an es- 
timate for r^d, M c i and p a v,ci with eqs. ([5}, <(Sj and (|10[) 
given the constraints on tm /t CIOSS , r^d/r^d and e found in 
the former section. 

Note that since q ow is a lower limit to the true SFE, 
and since the expansion factor, f r (res. gas), increases and 
the bound mass fraction, /m (res. gas), decreases with de- 
creasing SFE (figs. 1 and 4 of BK07), the estimated value 
for /, (res.gas) is an upper limit, and /jv/(res.gas) is a lower 
limit compared to their true values. In turn, our estimates on 
r h:Ci and M cl have to be seen as lower and upper estimates, 
respectively. 

4.3 Results 

Given the DMPP relation and a universal low-mass IMF, the 
results of our constraints (Sec. 14. 1} and the initial conditions 
for the DMPP GCs (Sec. S2J are as summarised in Tab. □ 

The estimated initial tidal-field strengths range from 
weak (rh,ci/r t ,ci ~ 0.02) to strong {r h ^i/r t>c i ~ 0.13). The 
left-most part of Fig.[3](c < 0.7) has only one observed clus- 
ter. This indicates that clusters forming in very strong tidal- 
fields may not su rvive for a Hubble-time an d are quickly 
tidally disrupted (jBaumgardt fc Makinoll2003h . 

We find a large range of gas expulsion time-scales 
for the observed clusters, from essentially instantaneous 
(TAf/icross <C 1) to adiabatic (r A f/t cr0 ss > 1) gas removal, 
which have to be explained by the constrained initial condi- 
tions (mass, size & metallicity). In Fig.[6]we plot the average 
metallicity for groups of clusters with similar gas removal 
time-scales (rjvi/t cross > 1,= 1 and< 1, respectively; com- 
pare Tab.Q]). We can thus confirm the idea from Sec.[2]that 
increasing metallicity reduces the time-scale of residual-gas 
expulsion and affects, in turn, the PDMF slope. 

Estimates to the lower limits of the SFEs range from 
flow ~ 10 to ~ 50 per cent, so essentially all SFEs used in 
the iV-body calculations occur. This covers observationally 
determined SFEs in young, g as embedded cluste rs which 
have SFEs of 20 - 40 per cent (jLada fc Ladall2003l l. 

2 Or two clusters, if either c \/r t c \ or tm Across matches the 
model grid. 



-1 



X 




relative residual-gas expulsion time x M /t cross 

Figure 6. The average metallicity for a group of clusters 
with similar constrained gas expulsion times is plotted against 
T M/*cross (Sec. [4731 . A clusters' metallicity determines whether 
the residual-gas is removed adiabatically (i.e. slow, Tjvf > tcross) 
or fast (rjv/ < tcross ) ■ Gas expulsion is more efficient in metal-rich 
systems, because the radiation driving the gas expulsion process 
couples more efficiently to metal-enriched gas, leading to larger 
expansion, mass-loss and, in turn, lower MF slopes (Sec. [2] Fig. 

rj. 



Most of the clusters share sub-pc initial half-mass radii 
(except Pal 5, NGC 288, NGC 6352, 6496 and 6712), sug- 
gesting that clusters form generally very compact, while Pal 
5 stands out with a very large birth radius (rh iC i = 7.52 pc, 
see discussion in Sec. I4.5|l . Initially small clusters have also 
today small half- mass radii (Tab. [2}. PD radii are ~ 3 — 10 
times larger than the birth values. Combining this with our 
constraints on the initial tidal-field strengths {rh,ta/i't,a\ m 
Tab. [TJ, we find the initial tidal radii, rt, to lie between 
~ 5 — 60 pc with an average of ~ 15 pc only. Such small 
tidal-radii imply that clusters appear to have formed un- 
der very extreme conditions. A possible qualitative picture 
of Galaxy formation, in which such environments may have 
been frequent, will be discussed in Sec. [5] 

Pre-cluster gas cloud-cores from which the stars of a 
cluster in our sample formed are estimated to have masses 
between - 3 x 10 5 M (NGC 6838) and ~ 2 x 10 7 M (NGC 
5139, i.e. ujCen). Initial cloud masses thus span a range of 
two orders of magnitude. Present day massive clusters gen- 
erally also had a massive progenitor cloud (Tab. [2}. Cloud 
masses were on average about 30 times larger than the PD 
cluster masses and about ~ 10 times larger than the Mfl n - 
values. The stellar mass, M ec i, was ~ 10 and ~ 2 times larger 
than the PD and M& n values, respectively. Pal 5, the cluster 
closest to dissolution (Sec. 14.5]) . had a ~ 50 x more massive 
progenitor cloud. NGC 5272s progenitor on the other hand 
was only ~ 12 x more massive than today. 

Initial average cloud densities range from - 10 2 (Pal 
5, has largest r ft , c .) up to ~ 10 8 M Q /pc 3 (NGC 7078, has 
very low ru,c\ and may have survived a SFE of 10 per cent). 
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Table 1. Present day and constrained initial values for the clusters in the DMPP sample. Columns denote from left to right: Catalogue 
number and cluster type (disc, young/old halo), low-mass MF slope, a, (modified) concentration, c and c mo( j feq. lllj l. tidal-field strength, 
r h,cl/ r t cli the time-scale of gas expulsion, tm /icross , the lower limit to the SFE, ei ow , the expansion factor and bound stellar mass 
fraction, f r = /, (res.gas) and /m = /M( res -g as )> half-mass radii, r h x , masses M x and average densities, p x = p aVl a;- The indices, x, 
deno te present day (PD), ini tial stellar (ecl=stars only) and cloud (cl=stars+gas) values. Afg n is the mass after residual-gas expulsion 
from lBaumgardt et all ibo08h . 

NGC type a c c mod r hjCl / r M I eiow fr fm r h>P T, r hjCl M PD M fin M ccl M cl p PD p e cl Pel 

[pc] [pc] [1O 5 M ] [1O 6 M ] [10 3 ] [10 6 M Q /pc 3 ] 

0.47 3.65 0.49 
0.39 5.36 1.63 
0.21 2.42 0.59 
0.45 19.98 7.52 
0.47 6.44 0.87 
0.40 3.38 0.47 
0.32 2.33 0.42 
0.29 3.07 0.64 
0.32 2.31 0.42 
0.53 2.60 0.23 
0.58 3.31 1.08 
0.61 1.55 0.17 
0.66 6.25 1.18 
0.61 3.03 0.42 
0.47 2.75 1.18 
0.53 2.72 0.24 
0.66 4.53 0.85 
0.56 1.87 0.47 
0.52 3.11 0.20 
0.53 2.67 0.23 
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Considering Pal 5 as an outlier, densities start at ~ 10 5 Mq 
pc -3 . Thus, densities decrease from cloud formation to their 
PD stellar densities by factors lying typically between ~ 10 2 
(NGC 6712) to ~ 10 5 (NGC 7078) with an average of ~ 
2 x 10 4 . The ratio of initial stellar to PD stellar densities is 
10 2 - 10 4 with an average at 3 x 10 3 . 



4.4 Combination of initial cluster parameters 

In this section we look for possible relations between the 
constrained initial cluster parameters. In order to see if and 
where a correlation between the parameters exist we com- 
piled Spearman rank-order correlation coefficients for the 
set of PD and initial cluster quantities used in this paper 
in Tab. [5] (upper half). We present or discuss those corre- 
lations which have the largest significant Spearman coeffi- 
cients (lower half of Tab. [2]). 

Young cl usters do not show any mas s-radius relation in 
observations (jLarsen 2004; Kroupa 2005) in agreement with 
anon-significant correlation between log(M cc i) and rh,c\- On 
the other hand their progenitor clouds are expected to show 
such a relation (Sec. [T]), but we see only a weak (maybe 
not significant) correlation between initial cloud half-mass 
radius and cluster mass from Fig. [7] However, we note the 
improvement of the correlation when comparing rh,d with 
Mci instead of A/ cc i which demonstrates the importance of 
comparing the relevant cluster parameters. 

We can understand the distribution of cluster size rec- 
ognizing its dependence on metallicity (Fig. [Sj . The initial 
half-mass radius appears to correlate with the cluster global 
metallicity in the sense that we find larger half-mass radii 
at larger metallicity. This trend can be expected, if higher- 
metallicity clusters can cool more efficiently and, thus, frag- 




0.1 1 10 100 

cloud mass M C | [10 6 M sun ] 

Figure 7. Mass of the cluster forming cloud, M c \, vs. initial half- 
mass radius, r h c i. Different cluster types are indicated. The clus- 
ter Pal 5 has c i = 7.52 pc at the mass indicated by the upward 
arrow. Most GCs share sub-pc initial half-mass radii. The half- 
mass radius appears to decrease with increasing cloud mass. The 
solid line is to guide the eye. 

ment earlier into individual stars. Metal-poor clusters have 
to wait for global collapse of the cluster to occur and a higher 
density before fragmentation becomes possible. Also, if the 
cluster formation process is feedback-regulated, as suggested 
by the existence of the maximum stellar-mass — star cluster- 
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Table 2. Spearman rank-order correlation coefficients, r s , between the quantities considered in the present paper (upper half). For no 
correlation, r s = 0, for perfect (anti-)correlation, r s = (— ) + 1. The probability, P(r s ), that the data show no correlation despite r s 
being different from zero, is given in the lower half of the table. We consider a relation to be significant, if P (r s ) 0.02. These numbers 
are bold-face in the lower part of this table. If numbers have a preceding '*', relations are considered to be trivial (e.g. mass and size 
are related to density, p oc M/r 3 ) or not of interest (why should the PD density be connected with the initial half-mass radius?). The 
correlations having a bold-face number in the upper half are depicted with a figure in this paper. The correlations in italic style are 
discussed. The rows and columns deno te: PD low-m ass MF slope; a; PD concentration, c; cluster metallicity, [Fe/H]; Galactocentric 
distance, i^gc (the latter two also from Harris 1996, 2003 revision); half-mass radii, r h x ; logarithmic values (lg = log 10 ) for the masses, 
M x , and average densities p x = pn V , x of the clusters from the DMPP sample. Indices as in Tab.fl] 



a c [Fe/H] d GC r hyPD r hzCl lg(A/ PD ) lg(A/ cc Q lg(A/ c i) lg(ppp) lg(Pcci) Ig(Pri) 

a 1.000 0.726 -0.647 0.104 - 0.069 - 0.733 0.646 0.191 0.598 0.109 0.783 0.804 

c 1.000 - 0.556 0.273 -0.405 - 0.848 0.540 0.367 0.680 0.450 0.931 0.911 

[Fe/H] 1.000 -0.365 0.247 0.636 -0.376 -0.164 -0.423 -0.295 -0.638 -0.657 

d GC 1.000 0.077 -0.133 -0.002 0.111 0.125 -0.024 0.192 0.167 

r h> p D 1.000 0.689 0.238 0.126 -0.003 -0.995 -0.505 -0.483 

r hyCl 1.000 -0.335 -0.131 -0.490 -0.719 -0.936 -0.930 

lg(M PD ) 1.000 0.741 0.870 -0.187 0.545 0.566 

lg(A4cl) 1-000 0.878 -0.078 0.397 0.402 

lg(M c i) 1.000 0.059 0.719 0.729 

lg(p PD ) 1.000 0.546 0.522 

lg(p C ci) 1.000 0.995 

lg(Pd) 1000 

a 0.000 0.000 0.002 0.663 0.771 0.000 0.002 0.421 0.005 0.649 0.000 0.000 

c 0.000 0.011 0.245 0.077 0.000 0.014 0.112 0.001 0.047 0.000 0.000 

[Fe/H] 0.000 0.113 0.295 0.003 0.102 0.490 0.063 0.207 0.002 0.002 

d GC 0.000 0.748 0.577 0.992 0.640 0.600 0.920 0.416 0.482 

r h>PD 0.000 0.001 0.313 0.596 0.990 *0.000 0.023 0.031 

r h>cl 0.000 0.149 0.581 0.028 *0.000 *0.000 *0.000 

lg(M PD ) 0.000 0.000 0.000 0.431 *0.013 *0.009 

IgGMeoi) 0.000 *0.000 0.743 0.083 0.079 

lg(M cl ) 0.000 0.806 *0.000 *0.000 

lg(pp D ) 0.000 *0.013 *0.018 

lg(pecl) 0.000 *0.000 

Ig(Pci) o-ooo 



mass relation (Wei dner. Kroupa fc Bonnelll l2009). then this 
correlation is naturally expected and would have a similar 
physical basis as metal-poor main-sequence stars being more 
compact than metal-rich stars. 

This view is supported when looking at initial cluster 
densities as a function of metallicity (Fig. [9}. Given our ex- 
planation for Fig. [8] we would expect and indeed see that the 
initial stellar (p av , C ci) and cloud density (p av ,ol) is larger in 
low-metalli city environments . The 2D hydrodynamical sim- 
ulations of iHocuk fc Spaanj (|2009h confirm the easier frag- 
mentation in more metal-rich environments and show that 
its dependence on metallicity is strong. 

Furthermore one expects gas expulsion to depend on 
the depth of the potential well the gas finds itself in (eq.Q]), 
which is confirmed by Fig. 1101 In a high-mass cloud the gas 
is removed less efficiently than from a low-mass one so the 
cluster retains more stars and a canonical a. Again, we point 
out that we see a correlation only because we looked at the 
progenitor cloud mass. The PDMF slope shows no relation 
with the total stellar mass forming out of the cloud material 
(Tab. El). 

The MF slope also correlates with metallicity. We have 
shown in Fig. [2] that the PDMF becomes more strongly 
depleted with increasing metallicity. Depleted and non- 
depleted observed clusters are distinctly apart in this di- 
agram. It is expected that high-metallicity gas is more eas- 
ily removed by radiation than metal metal-poor gas which 



ultimately leads to a PD depleted MF in initially mass- 
segregated clusters (Fig. |5J). 

Low-mass star depleted clusters have the largest initial 
half-mass radii (Tab.[2]| and, as a result, the lowest densities 
(Fig. 1 1 1 p . Large initial half-mass radii in low-a clusters will 
also lead to faster gas expulsion in them (eq. [1} . 

The PD concentration parameter, c, correlates with the 
same cluster parameters as a (Tab. which is obvious 
given the DMPP relation (eq. [2}. Thus, c also correlates 
strongly with the initial densities, even more strongly than 
with the PD densitjjf]. This suggests that the PD c-value is 
a measure for the initial density. 

4.5 Discussion 

The analysis above to constrain initial conditions of Galactic 
GCs yields remarkably nice results and information about 
how (globular) clusters might have formed. 

Based on the still small sample of clusters in Figs. [2] 
[8] and [9] metallicity appears to be an important parame- 
ter in determining the details of gas expulsion and cluster 
formation. The existence of metallicity trends with initial 
parameters further suggests that the observed metallicities 
are primordial and clusters can't be strongly self-enriched. 

3 Because c involves the radius ratio rt / r c it may not trace den- 
sity at all. 
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Figure 8. The initial half-mass radius, r h c i, increases with global 
cluster metallicity, [Fe/H]. The solid line presents a linear fit to 
all clusters in the sample. The dashed line is found excluding 
the disc and young halo clusters. The trend can be understood 
if metal-rich clusters cool more efficiently and fragment earlier 
into stars than metal-poor clusters, i.e. before global collapse of 
the star forming cloud occurs and also perhaps resulting from a 
feedback-regulated cluster formation process. Symbols as in Fig. 
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Figure 9. Initial average cloud density within the half-mass ra- 
dius, p 

av ,cl j SI 

function of cluster metallicity, [Fe/H]. Pal 5 has 
Pav cl ~ 10 2 Mq pc -3 at the metallicity indicated by the down- 
ward arrow. Metal-poor clusters are denser than the relatively 
more metal-rich clusters as expected from our explanation for 
Fig. [8] Symbols and best-fitting lines as in Fig. [8] 
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Figure 10. The PD global low-mass MF slope, a, is lower in 
star forming clouds with low masses, M c i, as expected from eq. [l] 
if residual-gas expulsion is the dominant source of low-mass star 
depiction. Symbols and best-fitting lines as in Figs. [8] and 
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Figure 11. Global low- mass MF slope, a, versus the density of 
the star forming cloud, p aVjC i. Pal 5 has p av c i ~ 10 2 Mq pc~ 3 
at the MF slope indicated by the leftward arrow. Low-mass star 
depleted clusters have the lowest concentration values today (Fig. 
[TJ and are also initially the least dense objects. Symbols and best- 
fitting lines as in Figs. 151 l9l and [TOl 



Metal-enrichment, however, occurs as seen by the multiple 
stellar populat ions observed in mas sive GCs. Next to NGC 
5139 (o;Cen. IVillanova et all 120071 ) and NGC 6656 (M22, 
iLee et al.|[2009l T7 Terzan 5 in the Galactic bulg e has recently 
been shown to host two stellar populations (|Ferraro et al.l 
2009). Metallicities for the different populations in uCen 
and Terzan 5 differ by A[Fe/H] = 0.5 dex at most. They 
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are both thought to be cores of former massive dwarf galax- 
ies such that they were able to retain their supern ova ejecta 
from which a second generation of stars formed l|Lee et al,l 
2009) and later merged with our Galaxy and can in that 
sense viewed to be exceptional. Alternatively, these clusters 
may have been immersed in a dwarf galaxy from which they 
repeatedly accreted in terstellar matter which initiated th e 
formation of new stars (|Pfiamm-Altenburg fc Kroupal2 009). 
In M22 the metallicity difference is much less, A[Fe/H] = 
0.2. Other clusters show virtually no dispersion in their iron 
content. That is, self-enrichment, if it did occur, did not 
significantly increase the bulk metallicity. 

We found that the primary parameter determining a 
is the initial density (the Spearman coefficient is largest, 
Fig. Hip . This results from a significant correlation of the 
PD a with the initial mass (Fig. I10[) and a significant anti- 
correlation with the initial size (Tab. [2]). The secondary pa- 
rameter influencing the PDMF slope is metallicity (Fig. [2]). 
As explained, both dependencies point to a link between 
a and the details of residual-gas loss. A dependence of a 
on the PD density does not exist and might indicate that 
two-body relaxation driven evolution is not as important as 
residual-gas expulsion in determining a. 

A lower initial density of the presently low-c clusters 
is also in concordance with the introduction of primordial 
unseen binaries in the analysis of MKB08: A binary frac- 
tion, /bin, was arbitrarily assigned according to the final 
c-value of the model clusters, i.e. after gas expulsion, with a 
larger value of /bin at low concentrations in order to better 
match the models with observations. This is now naturally 
explained since the initially denser clusters will have dis- 
solved more bi naries because interactions happen more often 
l|Kroupalll995h . The prediction of this work would thus be 
that presently low-c clusters ought to have a higher binary 
fraction. 

Neither initial nor PD parameters correlate with the PD 
Galactocentric distance, dGC, of the clusters (see Tab. [2}, 
although a weak anti-correlation between dec and [Fe/H] 
may be suggestive (Sec. [5]). While it is not clear why clus- 
ters should have for med at about the same distance where 
one sees them todav. lDiorgovskiT Piotto fc Capacciolil (Il993l ) 
found that the low-mass MF slope, a, is primarily deter- 
mined by dGC (and the height above the Galactic plane), 
indicating that the strength and change of the tidal-field 
(tidal shocking) in determining a dominates over internal 
processes. We can not confirm their results since we don't 
find a dependence of the PD a on dGC with the DMPP 
sample of clusters. The second pa rameter determin i ng the 
PD value of a in the analysis of iDiorgovski et all (l993) 
is metallicity, confirmed by our results. In their and subse- 
quent works this dependence was not understood beside the 
idea that the correlation has its origin at cluster formation. 
We already qualitatively excluded a varying low-mass part 
of the IMF in dependence of [Fe/H] to explain the DMPP 
relation since it is opposite to expectation (Sec. [2} and we 
instead suggest a metallicity dependent residual-gas expul- 
sion process to be the source of the depletion of low-mass 
stars in metal-rich, mass-segregated clusters (Figs. [2] and [6]). 

Pal 5 is outstanding in terms of its large initial ra- 
dius and low initial density. This might be the case because 
it is closest to dissolution, i.e. the effect of significant dy- 
namical evolution is important (see the next section). In 



this case the constraints might not be as good as for the 
other clusters. However, if the initial parameters for Pal 
5 were as constrained in this work, this clu ster might be 
the surviving remnant of a cluster complex (IKroupa 1998 



iFellhauer fc Kroup a 2002; B rims. Kroupa fc Felihauerll2009l 
see Sec. . 

For each of the figures presented in the former section 
we do not find significant differences between disc, young 
and old halo clusters. The best-fitting lines change some- 
what when excluding the disc and young halo clusters, but 
the improvement of the correlation measured by the Spear- 
man coefficients isn't strong. Despite the likely different con- 
ditions experienced by halo and disk GCs, respectively, at 
the time of their formation, especially due to different tidal- 
fields, this holds as evidence that the mechanisms are similar 
in the different types of GCs. As indicated by our constraints 
(Tab. [I] Fig. E}, the disc GCs NGC 6352, 6496 and 6838 
(except NGC 104) were born in stronger tidal-fields (larger 
rh,c\/r t ,c\) than most of the halo GCs. By the response of gas 
expulsion to the primordial tidal-field of their more extreme 
host environments, these clusters find their PD position in 
the DMPP diagram. 

Our constraints on the initial parameters rely, however, 
on the validity of the assumptions put into our analysis. 
The effects of dynamical and stellar evolution were included 
only in a simple way. A model taking gas expulsion and 
these effects self-consistently into account would lead to im- 
proved results: The dispersion seen in the figures including 
the PDMF slope, a, is relatively large (Figs. l2l [T0l and [Pi) 
and could possibly be significantly reduced. In the following 
we will therefore discuss the relative importance of dynam- 
ical and stellar evolution. 

4-5.1 Dynamical evolution 

The quality of the constraints on the initial tidal-field 
strength, the gas expulsion time-scale and the SFE (Sec. 
14 . 1 p depends on how strong the fingerprint of residual-gas 
expulsion on the PD values a and c is. Since MKB08 con- 
sidered only the effect of residual-gas expulsion but didn't 
include stellar and dynamical evolution, the change of the 
parameters over a Hubble-time doesn't enter their models. 
Secular evolution driven by two-body relaxation alters the 
low-mass MF slope by the preferential loss of low-mass stars 
and the concentration parameter due to processes such as 
core-collapse. We accounted for the latter by using the mod- 
ified concentration, c mo( j (eq. Ill|l . 

The ability of the iV-body models of BK07 to repro- 
duce the observations is a good starting-point to argue that 
residual-gas loss leaves a trace in the PD observables. This 
suggests that a and c are established rather quickly after 
residual-gas loss. Indeed, we have shown that the PD con- 
centration might as well be a direct result of the initial den- 
sity (Tab. [5J ■ And because of the dependence of PD a on 
M c i (Fig. [Toll, r h,d (Tab. [2]) and [Fe/H] (Fig. 0, respec- 
tively, we are also confident that residual-gas expulsion is 
an important step in initiating the observed PD value of a 
although this is the shortest, but on the other hand most 
violent phase in the life of a cluster; 

Further evidence is given by iBaumgardt et all (|2008h 
who explained the strongest low-mass star depleted clus- 
ters in the DMPP sample with secular evolution of mass- 
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segregated clusters without taking residual-gas expulsion 
into account. They showed that the observed PDMF slope is 
additionally a function of the remaining life-time of a cluster 
and calculated that the clusters that are only slightly defi- 
cient in low-mass stars (a ^ 1) still need more than or about 
a Hubble time until complete cluster dissolution. These clus- 
ters can be considered as dynamically young objects and our 
constraints are probably best for them. The low-mass star 
depleted clusters (a 1) have remaining lifetimes of 2 — 10 
Gyr with Pal 5 being closest to dissolution. 

Following our results, the low-a clusters are the least 
dense objects initially (Fig. II After residual-gas loss (and 
stellar mass loss, see below) they will have strongly fur- 
ther reduced their density. Indeed, MKB08 already consid- 
ered these objects as the fluffy remnants of residual-gas loss 
whose structural properties might not be altered strongly. 
So these clusters of the DMPP sample would possibly be 
least affected by dynamical evolution. 



4-5.2 Stellar evolution and mass-segregation 

As massive stars evolve they loose a fraction of their mass. 
When this mass is lost from the cluster, it will expand fur- 
ther changing its size an d enhancing the loss of stars over 
the tidal boundary (e.g. iBaumgardt fe Makind [20031). We 
included this effect by using the results of iHillsl (fl980l ) for 
the change of th e size (eq. [6]) and applying results from 
IBaumgardt et all J2008I) for the additional mass loss 



IS)). IVesperini. McMillan fc Portegies ZwartJ j2009l ) showed 



that the effect of stellar evolution is even larger in mass- 
segregated clusters and that such clusters, if they fill their 
tidal-radius, wou ld quickly disso l ve un der the influence of 
stellar mass loss. IVesperini et al.l j2009) conclude that clus- 
ters either can't have formed with a high degree of mass- 
segregation or they must be tidally-underfilled. Since gas 
expulsion in the models of BK07 naturally lead to tidally- 
filled clusters and MKB08 assumed strong mass-segregation, 
many of their model clusters would perhaps be quickly de- 
stroyed. Thus, the clusters observed today, as suggested by 
the results presented here, must have indeed formed highly 
compact and massive. 

To summarise, a lot of work remains to be done to fully 
understand the evolution of star clusters from their birth to 
dissolution. Therefore self-consistent numerical modelling, 
including all dissolution mechanisms, i.e. residual-gas expul- 
sion, stellar evolution and secular dynamical evolution, with 
realistic initial conditions (sizes and masses or equivalently, 
density, primordial binarity, ...) will be necessary, starting 
from the here derived initial conditions. 



5 ASSEMBLY OF THE GALACTIC OLD-GC 
HALO 

In this section we combine the previously found results and 
develop a picture of the formation of the population II, i.e. 
the metal-poor and old halo of GCs. 

We find that the gas expulsion time-scale, TM/t CIOSS , 
and the initial strength of the tidal-field, rh, c \/rt,ci, vary 
more strongly between low concentration clusters and, in 
turn, increase the scatter in the MF slope, a (Figs. [3] and 
[4]). Additionally, as the concentration, c, decreases, rh,ci/rt,ci 




Figure 12. Contraction scenario as described in Sec. [5] Left-hand 
picture: In a collapsing cloud the first clusters (the stars in the 
picture) form all over the cloud in a smooth potential. The condi- 
tions are similar for all of them leading to comparable MF slopes 
and concentrations after gas expulsion. These clusters enrich their 
immediate environment with metals. Right-hand picture: As the 
collapse proceeds more clusters form and the potential becomes 
grainy enhancing the mean tidal-field strength. Star clusters, e.g. 
the right filled cluster, forming near other clusters or dense clouds 
experience extreme conditions. They are more strongly enriched 
in metals from the surrounding objects than other clusters, e.g. 
the left filled cluster which is located in a more isolated place 
of the grainy potential, and they experience stronger tidal-fields 
than the first-forming clusters. This leads to strong differences be- 
tween the formation sites and, in turn, to variations in the PDMF 
slope. 



increases, i.e. PD low-c clusters have formed in stronger 
tidal-fields. The tidal-field strengths found for the DMPP 
clusters can't be explained with the PD configuration of the 
MW globular cluster system, i.e. a and the Galactocentric 
distance, dGC, show no correlation (see Tab. [2]). We argue 
that all these observations can be understood to be the re- 
sult of the assembly of the Galactic OH GC system within 
the initial cloud out of which the MW has formed. 

The OH clusters appear to be co-eval to a good ap- 
proximation dSalaris fc Weiss! |2002| ; |Pe Angeli et all 120051 ; 
M arm- Franch et al.ll2009l ). Despite their formation at the 
same time the clusters show a large spread in metallicity 
(-2.4 < [Fe/H] < -1.0 for the OH clusters in the DMPP 
sample). In Sec. l4.5l we argued that the metallicity in Galac- 
tic GCs might be primordial and has not been significantly 
altered. If the metal content of the more metal-rich old pop- 
ulation clusters comes from the products of stellar evolution 
which have been produced by earlier forming massive stars 
in other clusters, the inter-cluster medium (ICM) will be 
successively enriched with time. This material will then be 
recycled in later forming clusters. In this case a higher metal- 
licity as observed for the low-a GCs (Fig. [2} suggests that 
these clusters formed from an ICM enriched in metals and 
they should thus be somewhat younger than the relatively 
more metal-poor clusters. The age difference can then be at 
most a few hundred Myr only, corresponding to a dynamical 
time-scale. 

These ideas lead us to a picture of the formation of 
the old inner and co-eval Galactic GC system, which might 
have formed during the contraction of a giant gas cloud out 
of which finally, after an additional long period of ongo- 
ing accretion and me rging of formerly extragalactic systems 



ijSearle fc Zinnlll978h . the MW has emerged (possib ly along 
the lines originally proposed by lEggen et alj(l962l ). In situ 
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formation of the old GCs is in concordanc e with the age 
determination bv lMarin-Franch et alj (120091 1 who find that 
"the age dispersion of the old group of globular clusters is 
not in contradiction wi th the formation from the collapse of 
a single protosystem" . IZolotov et al.l (|2009h find the radial 
transition between an in-situ and accretion-dominated halo 
in MW like galaxies to lie around (fcc ~ 20 kpc, but most 
of the DMPP cluster even have rfG C < 10 kpc (e x cept NGC 
288, 2298, 5272, 7078 and Pal 5). ICarollo et~al1 (|2007i . and 
citations therein) argue that the dissipational merging of 
massive sub-galactic clumps builds up the inner halo which 
is followed by a stage of adiabatic flattening owing to the 
growth of the disk in order to explain the high-eccentricity 
orbits and higher peak metallicities observed for stars of the 
inner halo. While post-merger star formation may continue 
in the latte r phase, also GCs will for m in situ within the ini- 
tial galaxy. iForbes fc Bridges! i|2010r i point out that at least 
1/4 of the Galactic GC system are of extragalactic origin 
and that, although OH clusters are commonly thought to 
have formed in situ, some clusters with OH classifications 
might actually be accreted. 

However, whether still most the old halo GCs formed 
in situ, or whether the inner GC system is t he result of the 
merg ing of sub-galactic building blocks (e.g. Be kki fc Chibal 
2001), doesn't affect our discussion below. If the presented 
ideas are correct, we expect the scenario described there to 
be valid whereever the Galactic GCs formed, either during 
a single collapse or a collapse within the potentials of sub- 
structures which later merged, since the physical processes 
are fundamental. This can be supported by the two young 
halo clusters NGC 5272 and Pal 5, which possibly became 
part of the Galaxy during an accretion phase and fit per- 
fectly well into the DMPP diagram (Figs. 0© suggesting 
formation processes similar to those of old halo clusters and 
depending on the environment only (Sec. 14.51) . 

In this frame, clusters which are depleted of low-mass 
stars, i.e. those that reside in stronger tidal-fields (Fig. [3]) 
and expel their gas faster (Figs. [4] and [6]) because they are 
relatively more metal-rich (Fig. [SJ, were born at a later stage 
of (the GCs host-)galaxy formation. The environment in 
which the clusters formed must then have changed drasti- 
cally within a short time. 

The initial conditions can be understood, if the overall 
potential was rather smooth in the beginning. The first GCs 
formed all over the cloud experiencing similar smooth tidal- 
fields more or less independent of cJgc and explaining the 
comparable values of rh,ci/rt,ci- The global metal content 
was still low, so one finds generally moderate to slow gas ex- 
pulsion times (with respect to the crossing-time, Fig. [6]) and 
the clusters were smaller, i.e. more compact initially (Fig. [8]). 
These first clusters enriched the ICM of their local environ- 
ment with metals, from which the somewhat younger clus- 
ters of the old population were born. An increasing metal- 
licity with time leads to successively shorter gas expulsion 
time-scales (Fig. [6J and strong low-mass star loss in a mass- 
segregated cluster (Fig. [2]). The higher metallicity clusters 
were less compact than the previously formed clusters which 
had a lower metallicity (Fig. [9} and they suffered stronger 
tidal- fields than their metal-poorer counterparts (Figs.[2]and 

The pre-(GC-host-)galaxy gas cloud contracted due 
to self-gravitation during this process. The cloud may 



have become clumpy and substructures emerged. Frag- 
mentation of the cloud into massive star cluster form- 
ing regions made the potential grainy, thereby explaining 
the different and on average stronger tidal-field strengths 
among the younger clusters. Most of these clusters (~ 90 
per cent) will have been quickly dest royed by residual- 
gas expulsion (e.g. lLada fc Ladal 120031. and BK07) loos- 
ing their stars to the halo field jKroupa fc Boilvl |2002| ; 
iBaumgardt. Kroupa fc Parmentierl 12008! ) . Clusters forming 
next to other massive, dense objects are forced to expel their 
gas under more extreme conditions (a stronger tidal-field) 
than expected from their PD Galactocentric distance. 

Gas expulsion from early forming clusters may trig- 
ger the formation of new clusters (via gas compression) 
in their immediate surrounding perhaps leading to cluster 
complexes (CCs), i.e. clusters of star clusters. Such objects 
were recently found in the vicinity of spiral a rms in ex- 
ternal spiral galaxies (e.g. iBastian et al.1 120071 . and refer- 
ences therein). It is imaginable that similar processes were 
at work in the GC forming gas cloud. Observed CCs con- 
tain dozens to hundreds of st ar clusters and have radii up 
to a few hundred parsecs (e.g. IWhitmore fc~ Schwcizcr Il995l . 
for the Antennae galaxies). Merging processes in th e central 
dense parts of CCs may form more massive clusters dKroupal 
| l998 [; lFellhauer fc Kroupal2002l ; [Br"uns. Kroupa fc Fellhauerl 
2009). Clusters in such proximity may feel their mutual 
tidal-forces acting on each other and giving rise to small 
tidal- radii (as found for some DMPP clusters in Sec. 14. 3[) . 
Pal 5 may be such a surviving merged CC. 

All these effects lead to depleted MFs and lower con- 
centrations after gas expulsion in the younger of the OH 
clusters. In this sense the history of events that lead to the 
inner GC system may involve local rapid (on time-scales of 
hundreds of Myr) re-arrangements of the interstellar matter 
superposed on the overall contraction or collapse to the fi- 
nal population II spheroid. The weak anticorrelation between 
[Fe/H] and (fcc (Tab.[T}, though not very strongly signifi- 
cant according to our criterion, would be consistent with this 
scenario since enrichment continues while the cloud collapse 
proceeds and more-metal rich clusters would thus typically, 
but not exclusively, form closer to the centre. Our ideas of 
the formation of the Galactic GC system are summarised in 

Fig.m 



6 SUMMARY 

iMarks. Kroupa fc BaumgardtJ (|2008l , MKB08) showed that 
residu al-gas expulsion models by IBaumgardt fc Kroupa! 
(2007, BK07) explain the on first sight curious trend be- 
tween the global present day low-mass PDMF slope, a, in 
Galactic globular clusters (GCs) and their concen t ration 
parameter, c, (Fig. Q} as found by |Pe Marchi et al.l (120071 . 
DMPP): Weakly concentrated, low-mass-star-depleted clus- 
ters are a natural outcome of an initially mass-segregated 
cluster evolving through residual-gas expulsion. Given the 
DMPP data, it becomes possible to constrain the time-scale 
of gas expulsion, the tidal-field strength and the star forma- 
tion efficiency (SFE) all GCs in such a sample must have 
had. This is achieved by direct comparison of the MKB08 
models of residual-gas expulsion with the observations by 
DMPP and assuming the stellar low-mass IMF is universal 
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as is indicated to be the case from a large variety of stel- 
lar systems. The above parameters determine the expansion 
and mass loss following the residual-gas expulsion process 
in iV-body models. After applying corrections for the stellar 
and dynamical evolution of the clusters it becomes possible 
to calculate their initial half-mass radii, masses and densities 
directly before the residual-gas throw-out. 

We find that the primordial gas in different clusters is 
lost on short or long time-scales, measured with respect to 
the initial crossing-time of the cluster. Some GCs had SFEs 
as low as 10 per cent, others needed 50 per cent to sur- 
vive until today. Tidal-fields experienced by the clusters were 
weak or strong. Initial pre-cluster cloud-core masses are esti- 
mated to lie between 10 5 — 10 7 Mq and cluster radii typically 
were smaller than 1 pc. Average birth densities would thus 
have been 10 5 - 10 7 M pc~ 3 . 

We compiled Spearman rank-order correlation coeffi- 
cients to find connections between present day and initial 
cluster parameters and showed that: 

• Present day massive clusters also had massive progen- 
itors. Present day small-sized clusters initially also had a 
small half- mass radius (Tab. [2]). 

• The PDMF slope, the initial cluster size and density 
are functions of the PD cluster metallicity (Figs. [2] [8] and 

GUi. 

• Initial half-mass radii are larger and densities are lower 
in higher-metallicity environments, consistent with cooling 
being more efficient and fragmentation into stars occur- 
ring earlier before the overall global collapse of the cluster 
l|Hocuk fc Spaansll2009h . 

• The MF slope trend with metallicity (Fig. [2| could hint 
at a metal-dependent low-mass stellar IMF. However, the 
Jeans mass should be higher in low- metallicity environments 
and one would therefore expect a smaller fraction of low- 
mass stars there, contrary to observations, so that a varying 
low-mass part of the IMF seems not a viable explanation. 

• We interpret the PDMF observed to be depleted in 
low-mass stars in higher metallicity environments (Fig. [2]) 
as the result of a metallicity dependent gas expulsion pro- 
cess: Radiation from stars couples better to high-metallicity 
gas so that the residual-gas from star formation is expelled 
more efficiently and therefore probably quicker (Fig. [6]), sim- 
ilar to the me tallicity dependend mass-loss rates of evolv- 
ing stars (e.g. iMokiem et alJ [20071 ). This also indicates ef- 
ficient, i.e. fast, gas removal in clusters with a low value 
of the PDMF slope. If it is true, this provides for the first 
time an explanation for the dependence of the PDMF on 
the metallicit y, although it has been investigated in other 
works before jMcClure et al.lll986l : ISmith fc McClurdll987l : 
iDiorgovski et al.lll993l) .~ 

• Low-mass-star-depleted clusters have the lowest initial 
masses (Fig. I10|) and the largest initial radii (Tab. 0). This 
dependency of the PDMF slope may be due to more efficient 
gas throw-out in low-mass clusters with a large size (eq. 
[TJ. This then leads to stronger cluster expansion following 
residual-gas expulsion and, in case of a primordial mass- 
segregated cluster, to enhanced low-mass star loss across 
the tidal-radius. 

• As a result, the correlation between the density of the 
cluster-forming cloud-core and the PDMF slope, a, is found 
to be strongest of all (Fig. Ill[) . 



• We find no correlation of a with the distance of the 
cluster to the Galactic Centre (T ab. [2j , in contrast to the 
results of Djorgovski et all (|1993D . 

• The tendency of low-mass-star-depleted clusters to be 
initially least dense (Fig. lll|) allowed us to understand a sup- 
posed higher binary fraction in PD less-concentrated clus- 
ters: a higher fraction of unseen binaries had been arbitrarily 
introduced by MKB08 for weakly concentrated clusters to 
better match the observational data. Since the low-a clusters 
have today a low concentration (Fig. [TJ and are init ially less 
dense , fewer binary systems are dissolved in them (jKroupal 
Il995h and we predict a larger fraction of binaries in PD low-c 
clusters. 



Given the many correlations with metallicity, we pos- 
sibly unmasked metallicity as a main parameter in the for- 
mation process of clusters and for the residual-gas expulsion 
process. 

Knowledge about these initial conditions of star clus- 
ters lead us to the development of a picture of the for- 
mation of the old halo GC population within the frame- 
work of t he contraction of a giant pre-Milky W ay (MW) 
gas cloud (Egg en. Lvnden-Bell fc S andagc 1962). Since the 
clusters of the old halo are more or less co-eval and 
ages are con sistent with the formation from a single 



proto-system (ISalaris fc Weiss! 12001 : |Pe Angeli et al.ll2005l : 
M arin-Franch et al.ll2009l ). the age difference must be rather 
small (at most a few hundred Myr). From the initial con- 
ditions we see that these clusters have formed under very 
different conditions: Strongly concentrated, relatively metal- 
poor, and thus the oldest, clusters experienced weak tidal- 
fields only, i.e. the overall potential was smooth in the begin- 
ning so that conditions were rather similar in the pre-MW 
gas cloud. Gas was blown out of the clusters on compara- 
ble time-scales, because the metallicity was low (see above). 
Metal enrichment by these first clusters lead to successively 
higher [Fe/H] in later forming clusters. During that process 
the cloud contracted and became grainy on about a dynam- 
ical time-scale (few hundred Myr or less), perhaps triggered 
by residual-gas expulsion from other clusters and gas com- 
pression. 

Thus, by connecting the present-day concentration and 
low-mass stellar depletion of GCs it thus appears possible 
to uncover the evolution of the physical conditions during 
the first dynamical time of the MW: The environment 
in which clusters formed changed a lot as is seen in 
the stronger and also more strongly varying tidal-fields 
experienced by the weakly concentrated clusters that are 
typically more metal-rich. A larger global metallicity in 
clusters lead to more rapid gas blow-out and the easier 
loss of stars across the tidal-radius. This resulted in the 
younger, relatively metal-rich clusters appearing today as 
the low-concentration GCs with depleted low-mass stellar 
PDMFs. 
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